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1 Introduction 

Evolution equations of quantum many-particle systems arise in many problems of modern 
statistical mechanics. [1,10,12,18]. In the theory of such equations, during the last decade, 
many new results have been obtained, in particular concerning the fundamental problem 
of the rigorous derivation of quantum kinetic equations [19]- [27]. 

A description of quantum many-particle systems can be formulated in terms of two sets 
of objects: observables and states. The mean value defines a duality between observables 
and states and as a consequence there exist two approaches to the description of the 
evolution. Usually the evolution of many-particle systems is described, in the framework 
of the evolution of states, by the BBGKY hierarchy for marginal density operators. In 
the papers [21]- [27] a solution of the Cauchy problem to the quantum BBGKY hierarchy 
is constructed in the form of iteration series for initial data in the space of sequences of 
trace class operators. In [13, 14] for the quantum BBGKY hierarchy (for the classical 
many-particle systems in [15, 16]) a solution is represented in the form of series over 
particle clusters, whose evolution is described by the corresponding order cumulant (semi- 
invariant) of evolution operators of finitely many particles. Using an analog of Duhamel 
formulas, such a solution expansion reduces to an iteration series, which is valid for a 
particular class of initial data and interaction potentials. An equivalent approach for 
the description of many-particle system evolution is given by the evolution of observables 
and by the dual BBGKY hierarchy. For classical systems this approach is studied in the 
paper [2]. 

In this paper we consider the initial-value problem to the quantum dual BBGKY hier- 
archy of equations describing the evolution of observables for both finitely and infinitely 
many quantum particles, obeying Maxwell-Boltzmann statistics. We construct a solution 
of the Cauchy problem in the form of an expansion over clusters of the decreasing number 
of particles, whose evolution is governed by the corresponding-order dual cumulant (dual 
semi-invariant) of the evolution groups of operators of finitely many particles (groups of 
operators for the Heisenberg equations). 

We now outline the structure of the paper and the main results. 

In Section 2 we define the dual BBGKY hierarchy for quantum systems with a general 
type of an interaction potential between particles and construct a solution of the Cauchy 
problem to this hierarchy. We define the dual cumulants (the dual semi-invariants) of 
evolution operators for the Heisenberg equations and investigate some of their typical 
properties. The concept of cumulants of evolution operators forms the basis for the 
solution expansions of various evolution equations for infinitely many particles [13,16]. In 
Section 3 we study the properties of a one-parameter group of automorphisms generated 
by a solution of the initial-value problem to the dual BBGKY hierarchy in the space of 
sequences of bounded operators and prove the existence and uniqueness theorem of a 
classical and a generalized solution for the corresponding Cauchy problem. In Section 
4 we consider the description of evolution for infinite-particle systems in the Heisenberg 
picture and prove the existence of mean value functionals for observables. In the final 
section we give some concluding remarks. 

Let us introduce some preliminary definitions and results for the description of quantum 
many-particle systems of a non-fixed number of particles. 
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1.1 The initial- value problem to the Heisenberg equation 

We consider a quantum system of a non-fixed number of identical (spinless) particles with 
unit mass m = 1 in the space M?, v > 1 (in the terminology of statistical mechanics it is 
known as non-equilibrium grand canonical ensemble [1]). 

Let a sequence g = (I, gx, . . . , g n , . . . ) be an infinite sequence of self-adjoint bounded 
operators g n defined on the Fock space Tu = ©^Lo^^™ over ^ ne Hilbert space Tt (TC° = 
C and J is a unit operator). An operator g n , defined in the n-particle Hilbert space 
Tin — 7~C® n , will be denoted by g n (l, . . . , n). For a system of identical particles obeying 
Maxwell-Boltzmann statistics, one has g n (l, . . . , n) = g n {h, ■ ■ ■ , in) for any permutation 
of indexes {ii, . . . , i n } G {1, . . . , n}. 

Let the space £ i (J 7 h) be the space of sequences g = (I, g n , ...J of bounded 

operators g n (I is a unit operator) defined on the Hilbert space 7i n and satisfying symme- 
try property: g n (l, ...,n) = g n (ii, ...,i n ), if {h, ...,i n } G {l,...,n}, with an operator 
norm [5]. We will also consider a more general space £ 7 (JF W ) with a norm 

\\g\\^n) : =™ x ^ ll^lUc^), 

where < 7 < 1 and ||.|| £,m n ) is an operator norm [5]. 

An observable of finitely many quantum particles is a sequence of self-adjoint operators 
from Sj^T-h)- The case of the unbounded observables can be reduced to the case under 
consideration [3]. 

The evolution of observables A{t) = (I, A\(t, 1), . . . , A n (t, 1, . . . , n), . . . ) is described 
by the initial-value problem to the Heisenberg equation [3, 12] 

jA{t)=MA{tl (1) 
A(t)\ t=0 = A(0), (2) 

where, if g G V{H) C £(JF W ), the von Neumann operator Af = ®^l A/' n is defined by the 
formula 

(Ng)n ■= t {g n H n - H n g n ) . (3) 

h = 2irh is the Planck constant and the Hamiltonian H = ®^L -ff n in ([3]) is a self-adjoint 
operator with a domain of the definition T>{H) = {ip = @^ = Qip n G Tu I i>n G D{.H Tl ) C 

W n , E„ll^n|| 2 < OO} C ^ [5]. 

Assume Tt = L 2 (W), then an element ip G Tu = (B^qL 2 (R un ) is a sequence of func- 
tions ip = (ipo,^Pi (<?i), • • - ,i>n{qi, ■ ■ ■ ,Qn), ■ ■ •) such that ||^|| 2 = iV'ol 2 + J2n=l fdqi-. -dq n 
\ip n (qi, ■ ■ ■ ,q n )\ 2 < +00. On the subspace Ll(W n ) C L 2 (R un ) of infinitely different iable 
functions with compact support the n-particle Hamiltonian H n acts according to the 
formula (Hq = 0) 

j_2 n n n 

H n ij n = -~Y,\^n + Yl ZJ ® {k \qh,---,qi k )iJn, (4) 

i=l k=l «i<...<ifc=l 

where is a fc-body interaction potential satisfying Kato conditions [5]. 
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On the space £ / (J r n) a solution of the Cauchy problem to the Heisenberg equation (DO) 
is determined by the following one-parameter family of operators 

R 1 3t^g(t)g:=U(t)gU- l (t), (5) 

where U(t) = ®^L W n (i) and operators U n {t), U~ l {t) are defined as follows 

U n {t) :=ei tH ", 

U-\t) := e"**^ (6) 

and Uo(t) — I is a unit operator. 

Proposition 1. On the space £ 7 (JF W ) one-parameter mapping (J5]) 

M 1 3 t ^ g(t)g :=U{t)gU- 1 {t) 

defines an isometric *-weak continuous group of operators, i.e. it is a Cq -group. The 
infinitesimal generator M = ©^L M n of this group of operators is a closed operator for 
the *-weak topology and on its domain of the definition T>(N) C £ 7 (JF W ) which every- 
where dense for the *-weak topology, Af is defined as follows in the sense of the *-weak 
convergence of the space £ 7 (JF W ) 

w *-]hn- t (G(t)g - g) = -^(gH - Hg), (7) 

where H = @™ =0 H n is the Hamiltonian (j3j) and the operator: Mg = {—i/K){gH — Hg) 
is defined on the domain T>(H) C Tu- 

The validity of this statement follows from properties of one-parameter groups ([6]) [3] , 
[4]. Group of operators §5§ preserves the self-adjointness of operators. 

On the space £ 7 (^ 7 -^), for an abstract initial-value problem ([I])-©, the following state- 
ment holds 

Proposition 2. A solution of the initial-value problem to the Heisenberg equation (1T1)-(12"1) 

is determined by 

A{t) = g(t)A(0), (8) 

where the one-parameter family {Q(t)} t eR of operators is defined by expression (JSj). For 
A(0) E V(Af) C £ 7 (JF W ), A(t) is a classical solution and for arbitrary initial data A(0) G 
£ 7 (jF w ), A(t) is a generalized solution. 

According general theorems about properties of the dual semigroup [4], [6], [7] the 
validity of these propositions follows also from the fact that the group {{?(£) }<gR is dual 
to the strong continuous group of the von Neumann equation defined on the space of 
sequences of trace class operators [3]. 

We remark that, in the framework of kernels and symbols of operators, infinitesimal 
generator ([7j) and groups (jSJ) are studied in [11]. 
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1.2 Quantum many-particle systems 

A positive continuous linear functional on the space of observables A(t) E £(.F W ), whose 
value is interpreted as its mean value (average value of an observable), is defined as 
follows [4] 

OO ^ OO j 

(A(t)\D(0)) := ( ]T -Tr lv .., n D n (0)) ~ X ]T -Tr lr .., n A n (t, 1, . . . , n) D„(0, 1, . . . , n), (9) 



n! 14 — ' ni 

n=0 n=0 



where Tri r . , jn is the partial trace over 1, . . . , n particles and Yu^o^i^h-,^ D n (0) is a nor- 
malizing factor (grand canonical partition function). The sequences D(0) = (l,Di(0, 1), 
. . . , D n (0, 1, . . . , n), . . . ) of positive self-adjoint operators D n , n > 1, (the density op- 
erators whose kernels are known as the density matrices [11]) defined on the n-particle 
Hilbert space Ti n = TL® n = L 2 (R un ), describe the states of a quantum system of non-fixed 
number of particles. 

Usually it is assumed that the states belong to the space (.F^) = @^1 ot n £}{l-L n ) 
of sequences / = (I, f\, . . . , f n , . . .) of trace class operators /„ = / n (l, . . . , n) E £ 1 (H n ), 
satisfying the Maxwell-Boltzmann statistics symmetry condition, equipped with the trace 
norm 

OO OO 

H/lk(^) : =E «1/»II*CH») : =E « n Tr w |/ n (l,...,n)|, 

n=0 ra=0 

where a > 1 is a real number, Tri n is the partial trace over 1, . . . , n particles. We 
will denote by the everywhere dense set in £^,(JF W ) of finite sequences of degenerate 
operators (finite-rank operators) [5] with infinitely differentiate kernels with compact 
support. The space £^(JF W ) contains sequences of operators more general than those 
determining the states of systems. We will also consider the space £ 1 (JF W ) = ©^ ^ 1 (^)- 

If A(t) G £ 7 (.F W ) and D(0) E S^FjA, then functional © exists. 

It is well known that there exists an equivalent approach to the description of the 
evolution of quantum systems, which is given by the evolution of states. In this case the 
evolution of all possible states D(t) = (/, D\(t, 1), . . . , D n (t, 1, . . . , n), . . . ) is described by 
the initial- value problem to the von Neumann equation [3, 12] 

jD(t) = —AfD(t), (10) 
D(t)\ t=0 = D(0), (11) 

where if / G £q(JF w ) c V(J\f) C £ 1 (^ r -^) the von Neumann operator (—A/") is defined by 

{-Mf) n := ~{H n f n - f n H n ) , (12) 

The von Neumann equation (JTUl) is dual to the Heisenberg equation ([1]) with respect to 
bilinear form ([9]). 

For initial data -D(O) from the space £ 1 (JF W ) the solution 



D(t) = G(-t)D(0) 
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of the Cauchy problem to the von Neumann equation (TTO l) - ( TTTT) is determined by the 
following one-parameter family of operators on £ 1 (jF w ) 

R 1 Bt»g(-t)f:=U(-t)fU- 1 (-t), (13) 

where the operators U n (—t), IA~ X [— t) are defined by (jSJ). 

The properties of a one-parameter family {Q(—t)} t eR of operators (fl3"|) follow from the 
properties of groups (J6]). 

Proposition 3. On the space one-parameter mapping (|13p 

l^tM G(-t)f := U(-t)fU-\-t) 

defines an isometric strongly continuous group of operators, i.e. it is a Co-group, which 
preserves positivity and self-adjointness of operators. If f e £q(JF w ) c V(—Af) C £ 1 (JF W ) 
then the infinitesimal generator: —M = @^1 (— M n ) of this group of operators is deter- 
mined in the sense of the norm convergence in the space £ 1 (JF H ) as follows 

\^(g(-t)f-f) = - l -{Hf-fH) = -Nf, (14) 

where H = @^ =0 H n is the Hamiltonian (@J and the operator: (—i/h)(Hf — fH) is 
defined on the domain T>(H) C Tu- 

Thus, if we consider group ([5]) as the dual to group (fTBl . from the above Proposition 
it follows the statement about properties of group ([5]). 

For a system of a finite average number of particles there exists an equivalent possibility 
to describe observables and states, namely, by sequences of marginal observables (the so- 
called s-particle observables) G(t) = (G , Gi(t, 1), . . . , G s (t, 1, . . . , s), . . . ) and marginal 
states (or s-particle density operators) F(0) = (/, Fi(0, 1), . . . , -F s (0, 1, . . . , n), . . .) [1], 
[12]. These sequences are correspondingly introduced instead of sequences A(t) and D(0), 
in such way that mean value ([9]) does not change, i.e. 

oo 1 

(A(t)\D(0)) = (G(t)\F(0)) = - Tr i,..,s G B (t, 1, . . . , s)F a (0, l,...,s). (15) 

8 = S - 

Then marginal observables are defined by 

G s (t,Y):=J2-^f- E A s „ n (t,Y\{ 3l ,..., Jn }), s>l, (16) 

n=0 ' j 17 t...^j n =l 

where Y = (l,...,s) and the sequence A(t) = (i, Ai(t, 1), . . . , A n (t, 1, . . . , n), . . . ) is 
solution ([8]) of the initial- value problem to the Heisenberg equation ([I])-©, i.e. A{t) = 

g(t)A(o). ' 

If we introduce the operator a + (an analog of the creation operator [2]) 

{a+g) s (Y) := £ g^(y\{j}) (17) 



Quantum Dual BBGKY Hierarchy 



8 



defined on £ 7 (.F W ), then expression ( fl6l) can be rewritten in the following compact form 

G(t) =e- a+ g(t)A(0). 

The evolution of marginal observables (fTBj) of both finitely and infinitely many quan- 
tum particles is described by the initial- value problem to the dual BBGKY hierarchy. For 
finitely many particles the quantum dual BBGKY hierarchy is equivalent to the Heisen- 
berg equation flTJ) (the dual equation to the Heisenberg equation (pQ) is the von Neumann 
equation (TlOT) ) . For systems of classical particles the dual BBGKY hierarchy was examined 
in [1], [2], [15]. 



2 The quantum dual BBGKY hierarchy 

In order to approach the description of observables of quantum many-particle systems by 
the marginal observables (s-particle observables) we study the evolution of the system 
by means of the quantum dual BBGKY hierarchy. We introduce such a hierarchy of 
evolution equations and construct a solution of the Cauchy problem to this hierarchy. 



2.1 The initial-value problem to the quantum dual BBGKY 
hierarchy 

The evolution of marginal observables is described by the initial-value problem for the 
following hierarchy of evolution equations 

^-G s (t,Y)=M s (Y)G s (t,Y)+ (18) 
at 

E^ E jkh^ E ^Kh,...,j h )G s - n {t,Y\{j l ,...,j n }), 



G s (t) \ t=Q = G s {0), s > 1, (19) 
where, on P(A4) C £(Hk), the operator M^} is defined by 

^g n -.= ~{gn& n) -& n) g n ). (20) 



and the operator is introduced in Hamiltonian (jl]). We refer to f|T8|) as the quantum 
dual BBGKY hierarchy since the canonical BBGKY hierarchy [1] for marginal density 
operators is the dual hierarchy of evolution equations with respect to bilinear form ( TT5T) 
to evolution equations (fT8l) . 

In the case of two-body interaction potential (jl]), hierarchy (|T8|) has the form 

±G s (t,Y)=Af s (Y)G s (t,Y)+ tf(jij 2 )G s - 1 (i,y\{ji}), s > 1, (21) 
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where the operator J\f^ is defined by ( 120]) for n = 2. For H = L 2 (R^), the evolution of 
kernels of operators G s (t), s > 1, for equations (12T]) . is given by 

, q 3 ; q[, ■ ■ ■ , q' s ) = ( - j X)(-A ft + A ? ,) + 

gj) - $ (2) fe - qj)))G s {t, q 1 ,...,q s ;q[,...,q' s ) + 

q'j) - $ (2) (g, - qj ))G s ^(t, q u . . . , V, . . . , q s ; q[, . . . , V, . . . , q' s ), 

where (gi, . . . , V, . . . , q s ) = (q%, . . . , qj-i, qj+i, ■ ■ ■ , q s )- The dual BBGKY hierarchy for a 
system of classical particles stated in [1,2] is defined by similar recurrence evolution 
equations. 

The quantum dual BBGKY hierarchy ffTHj) can be derived from the sequence of the 
Heisenberg equations ([I]) provided that observables of a system are described in terms of 
marginal operators (s-particle observables) (|T6|) . 

Remark 1. Another way of looking to the derivation of the quantum dual BBGKY hier- 
archy consists in the construction of adjoint (dual) equations with respect to the bilinear 
form (TT5]) to the quantum BBGKY hierarchy. For the sequence F(t) = (/, Fi(t, 1), . . . , 
F s (t, Y), . . . ) of s-particle density operators (marginal density operators) F s (t,Y), s > 1, 
the quantum BBGKY hierarchy has the form [13] 

j/s(t,Y) = -Af s (Y)F s (t,Y) + 

S _ S CO - 

E^! E E^«.-^(-^£ +n) ) ,ik,X\Y)F a+n (t,X), (22) 

fc=l ' i 1 ^z...z/zi h =\ n=l 

where X = {1, . . . , s + n} and on £,l(H s+n ) C 2}(J-L s +n) the operator A/" i ^ +?1 ' ) is defined by 
(|20|) . Indeed, hierarchy (TT8I) is dual to hierarchy (|22|) with respect to bilinear form (|T5|) . 

Remark 2. In the paper [2], for classical systems of particles with a two-body interaction 
potential, an equivalent representation for the dual hierarchy generator was used. In 
the case under consideration, on the subspace £>(Q3 + ) C £ 7 (^-^), the generator has the 
following representation 

23+ = Af + [AT, a + ] , 

where [. , .1 is a commutator and the operator o + is defined on the space £ 7 (JF W ) by 
(1T7|) . In a general case the generator of the quantum dual BBGKY hierarchy (1181) can be 
represented in the following form 

23 + =A/>]T-i[...[A/;a + ],...,a + ] =e~ a+ Ne a \ (23) 

n-times 



d 

iU—G s {t,q u ... 



E ( $(2) ^' 



i=i<j 



E ( $(2) (^ 
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where the operators e are defined on the space £ 7 (JF W ) by the expansions 

(^).M = E^T ^ 9s-n(Y\{jx,...Jn}), 8>1. 

Representation (1251) is correct in consequence of definition (|T7|) of the operator a + and 
the validity of identity 

([...[AT,a+]__a^],) s (r) = 

n-times 

S j s 

E E ^i2 ) (jl 5 --- 5 Jfc)^-n( ir \Ol 5 ---:Jn}), 

fc=n+l ^ r j l9 t...^ij k =l 

which for a two-body interaction potential reduces to the following one 

s 

{[Af,a + ]g),(Y) = ^St(ji,h)g s -x{Y\{j 1 }). 

2.2 The formula for a solution and its representations 

We consider two different approaches to the construction of a solution of the quantum 
dual BBGKY hierarchy ( [18]) . Since hierarchy ( Tl8l) has the structure of recurrence equa- 
tions, we deduce that the solution can be constructed by successive integration of the 
inhomogeneous Heisenberg equations. Indeed, for solutions of the first two equations we 
obtain 

G 2 (t,l,2) = 2 (t,l,2)G 2 (O,l,2) + (24) 



dhg 2 {t -h,l, 2yv2(l, 2)0i (ti, l)0i (ti, 2) (Gi(0, 1) + Gi(0, 2)) . 

Let us transform the second term on the right hand side of (1241) for G 2 (£) as follows 

t 

J dt x Q 2 (t - tx, 1, 2).A/2(1, 2)0i (ti, l)0i (tx, 2) (25) 
o 

t 

/d 
dtijfi&it ~ h, 1, 2)0i (ti, l)0i {h, 2)) = 2 (t, 1, 2) - 0i (t, l)0i (t, 2). 

o 

The operator 2 (t, 1,2) - 0i(t, l)0i(£, 2) := 2l+(t,l,2) in equality (j2S} is the 2nd-order 
dual cumulant of evolution operators £?(£) = ©^l o 0n(£) ©• Equality (1251) is an analog 
of the Duhamel formula [9], which holds rigorously, for example, for bounded interaction 
potentials. 

Thus, for solutions of the first two equations of hierarchy ( fl8l) we finally obtain 
Gi(t,l) = 2l+(t,l)Gi(0,l), 

G 2 (t, 1, 2) = 21+ (t, ({1, 2})i)G 2 (0, 1,2)+ 21+ (t, 1, 2) (Gi(0, 1) + Gi(0, 2)) , 
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where we introduced the lst-order dual cumulant %lf(t, ({1, 2})i) := Qiit, 1,2) and the 
notation ({l,2})i denotes a set composed by only one of elements 1 and 2 and it is 
similarly generalized in next formula (j26|) . 

Making use of transformations similar to (125|) . for n > 2, a solution of equations 
(|T8|) . constructed by successive integration of the inhomogeneous Heisenberg equations, 
is represented by the following expansions 

s 1 s 

G s (t,Y) = J2j 7 z^_ E %t +n {t,(Y\X) 1 ,X)G s _ n (0,Y\X), s>l. (26) 

Here we used some abridged notations: K = (1, . . . , s), X = Y\{ji, . . . ,j s _ n }, the set 
(Y\X)i consists of one element from Y\X = (ji, ■ ■ ■ ,j s -n), i-e. the set (ji, . . . ,j s - n ) is 
a connected subset of the partition P (|P| = 1, |P| denotes the number of considered 
partitions). The evolution operator %lf +n (t, (Y\X)i, X) in fl26|) is the (1 + n)th-order 
dual cumulant 2l^ +n (t), n > 0, of groups of operators (J5]) 

%t +n {t,(Y\X)i,X)--= E (-1) |PM (|P|-1)! n^l(*' X «)' ^ 

p : {(y\^)i.x}=U^ XiCP 

where J2p * s the sum over all possible partitions P of the set {(Y\X) 1 , ji, . . . ,j s - n } into 
|P| nonempty mutually disjoint subsets Xi C {(Y\X)i, X}. We consider some properties 
and examples of dual cumulants fl27j) in the next subsection. 
Using the identity 

S ^ S S j s 

E ( s - n )\ E 9s-n(h,---,j.-n)=^2— l E ^-n(y\0'l,--.,Jn}) (28) 

which is valid for the Maxwell-Boltzmann statistics symmetry property, expansion ( j26l) 
can be rewritten in the following an equivalent form 

G.{t,Y) = 
s -i s 

E ~\ E *tn(t, {Y\{Jl, ■ ■ ■ M) V jl, ■ ■ ■ ,3n) G s -n(0, Y\{j,, . . . ,j n }), S > 1. 
n =0 ' ji^...^j„=l 

The formula for a solution of the quantum dual BBGKY hierarchy f[T8"j) can be also 
derived from solution ([8]) of the initial- value problem to the Heisenberg equation ([I])-© 
on the basis of expansions (TLB"]) . 

For sake of simplicity we consider the additive-type observables in capacity of initial 
data, i.e. one-component sequences G^(0) = (0, cti(l), 0, . . . ). In this case the expansion 
for solution f[2"61) attains the form 

s 

GV(t,Y)=Vlf(t,l,...,s)Y,*iti)i s>l. (29) 

To determine the unknown evolution operator 21+ (t) from expansion (1291) . we equate 
expansion (1291) for the operator G^ (t, Y) with its representation by formula (TIB"]) for the 
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additive-type observables A(0) = (0, ai(l), ■ ■ ■ ,Y1 a i(0> ■ ■ ■)> i- e - 

a+(t,i,...,s)^ ai (j)= (30) 

™=o ' ii^-^in=i iey\{ii,...,i„} 

Solving recursion relations (1301) for arbitrary operators X^=i Q iO) we fi 11 ^ that the 
order dual cumulant 2l + (t, 1, . . . , s) is defined by expansion (I34I) . which is particular case 
of (1271) . For example, the first two equations from recursion relations (1301) have the form 

2tf (M)ai(l) = &(M)ai(i), 

21+ (f, 1, 2) (oi(l) + oi(2)) = G 2 (t, 1, 2) ( 0l (l) + ai(2)) - 1)^(1) - 2)a x (2) = 
(&(t, 1, 2) - &(t, l)&(t, 2)) (ai(l) + oi(2)). 

Hence a solution of the quantum dual BBGKY hierarchy (1181) is defined by expansion 
( 126]) where the evolution operators 2l + +n (t), n > 0, are cumulants (j2"Tj) of groups (jSJ). 

Remark 3. In the paper [2], for classical systems of particles, an equivalent representation 
for the dual hierarchy solution was used. For the case under consideration the solution 
expansion has the following representation 

oo 1 

G(t) = Q(t)G(0) + [G(tl * + l---,* + ] G(0) = (31) 

n-times 

e- a+ G(t)e a+ G(0), 

where [. , . 1 is a commutator, the operator a + is defined on the space £ 7 (.F W ) by (ITTj) 
and the one-parameter family {<?(£) }t e R of operators is defined by expression (J5j). 

Representation (131!) is true in consequence of definition (fTTj) of the operator a + and 
the validity of identity ( 1281) and of the following 

(^[...[g(t),a^^]g) s (Y)= (32) 
n-times 

— ^- E X£ ¥n (t,{Y\X) 1 ,X)g 4 ^Qr\X), s>l, 
where X = Y\{j\, . . . ,j s - n } and Y\X = (ji, . . . ,j s - n )- For example, if n = 1 we have 

s 

([£(t),a+]s) s (r) = £ (&(t,y) -a s -i(*,nO'}))^-i(nO'}) = 

3=1 

j2^t^( Y \{j})i,j)9s-i(Y\{j}). 

3=1 

Some details of the calculations for identity (|32|) will be given in a general case in the 
Remark 4 of next subsection. 



Quantum Dual BBGKY Hierarchy 



13 



2.3 Dual cumulants of groups of evolution operators 

We now consider some properties of dual cumulants (|27|) . Let us expand the groups of 
operators Q{t) = ©^L {?n(i) © over new evolution operators as the following cluster 
expansions 

Gn(t,Y)= £ R %+(t,Xi), n=\Y\>0, (33) 
p : y=U ! x l x,cp 

where ^ p is the sum over all possible partitions P of the set Y = (l,...,n) into |P| 
nonempty mutually disjoint subsets X, C Y. A solution of recurrence relations ([3"31 is 
determined by the expansions [15] 

K(f,Y)= J2 (-l) |PM (|Pi-l)! II SmfrXi), n=m>0, (34) 

where the notations are similar to that in ([3"3"j) . We refer to the evolution operator 2l + (t) 
as the nt/i-order dual cumulant (dual semi- invariant) of group of operators (jSJ). 
The simplest examples of dual cumulants (1341 have the form 

2l+(M) = £i(M), 

a+(t, i,2) = g 2 (t,i,2)-g 1 {t, i)Gi{t,2), 

21+ (t, 1, 2, 3) = ^ 3 (t, 1, 2, 3) - &(t, 3)&(t, 1, 2) - 2)&(t, 1, 3) 
-Gi(t, l)g 2 (t, 2, 3) + 2!&(t, 1)& (t, 2)ft(t, 3). 

In the case of a quantum system of non-interacting particles (n > 2) we have: 2l + (t) = 
0. Indeed, for non-interacting quantum particles obeying Maxwell-Boltzmann statistics it 
holds 

n 

Q n (t, l,...,n)= JJ</i(t,i) 

i=i 

and hence 

Kit,Y)= y: (-i) |phi (ipi-i)'nnr>^^) 

P:Y=\J t X t X t CP 
n n 

= J2(- 1 ) k ~ 1 <n,k)(k-l)\l[g 1 (t,i) = 0, 

k=l i=l 

where s(n, fc) are the Stirling numbers of the second kind. Here the following equality is 
used 

n 

E (- 1 ) |Phl (i p i - 1 ) ! = E(- 1 ) fe_ls ^> - x ) ! = Vi. ( 35 ) 

P:Y=\J.Xi fc=l 

where 5 n> i is a Kroneker symbol. 

The generator of the lst-order dual cumulant is given by operator ([3]), i.e. 

w*-limi(2l + (t,F) -l)g n (Y) = Af n (Y)g n (Y), 

where for g n G T>(J\f n ) C £(7Y n ) this limit exists in the sense of the *-weak convergence 
of the space £(H n ) (Proposition 1). 
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In the general case an infinitesimal generator of the nth-order dual cumulant, n > 2, 
is the operator (A^ ) defined by n-body interaction potential (jl]). Indeed, according to 
equality (l35j) in the sense of a point-by-point convergence in the space £(TC n ) for the 
nth-order dual cumulant, n > 2, we have 

hm- t K(t,Y)g n (Y) = E ( — 1 ) 1 P 1 " 1 ( | P | - 1)!^ = 

p : y=U ! ^ ! XiCP 

E (-i) |pm (|pi-i)'E E E (,!,..., ^)^(r), 

P-.Y^Xi XiCP k=2 i 1< ...<i k e{X l } 

where Aff^ is defined as in ( 1201) . i.e. 

where the operator $( n ) is defined in Hamiltonian (T4j). 

Summing the coefficients of every operator Af^t we deduce 

w*-\imfy(t,Y)g n (Y) = M&\Y)g n (Y), (36) 

Thus, if g n G V(Af n ) C £(7Y n ), the generator of the nth-order dual cumulant is defined 
by (1361) in the sense of the *-weak convergence of the space £(H n ). 

We remark that, at the initial time t = 0, solution (I26p satisfies initial condition (1191) . 
Indeed, according to definitions ([HD (i^ =1 (0) = / is a unit operator) and of equality (1551) 
for n > 1, we have 

^ + n(0, (Y\X) h X) = E ("l) |Phl (|P| - W = ^n,l- 

P:{(Y\X)i,X}=U i ^ 

Remark 4. Cluster expansions (133]) can be put at the basis of all possible solution 
representations of the quantum dual BBGKY hierarchy (1T8|) . In fact we can obtain 
representation (J3T1) solving recurrence relations (J33l) with respect to the lst-order dual 
cumulants for the separation terms, which are independent from the variable Y\X = 
(ji, . . . , js- n ) 

Qltn&iYXX^X) = J2%t{tiY\XUZ) E (- 1 )' P| |P|! lC*t{t,Xi)> 

ZCX P:X\Z=|J l X 1 

where J2zcx * s a sum over a ^ subsets Z C X of the set X. Then, taking into account 
the identity 

E (- 1 ) 1151 l P l' lt2^^ X i)9s~n(Y\X) = E (- 1 ) 1151 iPN-n^) 

and the equality 

E (-l) |P| |P|! = (-l) |X ^, (37) 
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for expansion ( l26i) we have 

S -. s 

G.(t, Y) = J2 773^7 E E &-*4-[z| (*, Y\X U Z) G s _ n (0, Y\X). 

Thus, as a result of the symmetry property for the Maxwell-Boltzmann statistics and of 
definition (ITT)) of the operator a + , we derive representation (I3T1) 

e- a+ G(t)e a+ G(0). 

We can obtain one more representation for a solution of the initial-value problem to 
the quantum dual BBGKY hierarchy if we express the dual cumulants 2li+ n (t), n > 1, 
of group of operators (jSJ) with respect to the lst-order and 2n<i-order dual cumulants. In 
fact it holds 

%l +n {t,(Y\X) u X) = E *+(t,Y\X,Z) E (" 1 ) |P| l p l ! n!l'^(*'^)> 

P:X\Z=\J t X t 

Z^% 

where Ylzcx, is a sum over all nonempty subsets Z C X of the set X. Then taking into 
account the identity 

E (" 1 )' P| l P l' lt^^ X i)9s-n{Y\X) = E (" 1 ) |P| n9s-n(Y\X) 

and equality (1571) . we get the following representation for expansion (1251) 

G,(t,y) = s^(t,y)G,(o,y) + 

E E E ^^^^^(O,^), 

n=l V >' j 1 ^...^j s _ n =l ZCX, 

Z^V) 

where Y = (1, . . . , s), X = Y\{ji, . . . , j s _„}, i.e. Y\X = (j u . . . , j a _ n ). 

3 A solution of the initial- value problem to the 
quantum dual BBGKY hierarchy 

One-parameter mapping generated by solution ( l26l) is not strong continuous on the space 
£"y(J~n). Thus, the constructed solution satisfies the initial- value problem to the dual 
BBGKY hierarchy (TT8l)-(|T9l) only in the sense of the *-weak convergence of the space 
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3.1 A group of operators for the quantum dual BBGKY hier- 
archy 

On the space £ 7 (jF w ) solution (1261) of the initial- value problem to the dual BBGKY hier- 
archy (jl8j) - (ll9j) is determined by a one-parameter mapping with the following properties. 

Theorem 1. If g G £ 7 (^-^) and 7 < e~ x , then the one-parameter mapping 

(U + (t)g) s (Y) := (38) 

S j s 

EtJ^T E *l +n (t,(Y\X) 1 ,X)g s _ n (Y\X), s>l 

is a Cq -group. The infinitesimal generator QS + = ®^1 23^ of this group of operators is 
a closed operator for the *-weak topology and on the domain of the definition T>( ( B + ) C 
£ 7 (^ r - w ) which is the everywhere dense set for the *-weak topology of the space £ 7 (J^) it 
is defined by the operator 

(<B + g) s (Y):=Af s (Y)g s (Y) + (39) 
1 1 

E^i E (k-n)\ £ ^t(h,---,jk)9s-n(Y\{j 1 ,...,j n }), 

n=l • k=n+l 1 '■ j 19 t...ft k =l 

where the operator Af^ is given by fTSOj) . 



Proof. If g G £ 7 (JF W ), mapping (1381) is defined provided that 7 < e 1 and that the 
following estimate holds 

II^MUc^) ^ e2 ( X - T^^lblkc^)- ( 4 °) 
This estimate comes out from the inequality 

IK^^IU) < E E E (i p i - i)'ii^-niu (Ws _„) < 

n=0 V >' jijt...ft t - n =l P:{y\X)i,X}=|J i ^< 

s I n+1 

E ll^«-n||fl(W.-n)-f7~~ ^7 E S ( n + X » ~ 

n=0 v ' k=l 

where s(n + 1, fc) are the Stirling numbers of the second kind for which it holds 

n+1 n+1 ^ / , «\i n+1 

£.(„ +!,*)<* -!>! = £ I E ±i±^<yj*»< n!e »« 

fc=l fc=l ri,...,rfe>l " fc=l 

ri+...+r fe =n+l 

On the space £ 7 (JF W ) the *-weak continuity property of the group U + (t) over the 
parameter t G M 1 is a consequence of the *-weak continuity of group of operators <^ for 
the Heisenberg equation ([!]). 

In order to construct an infinitesimal generator of the group {U + (t)} t€WL we firstly 
differentiate the nth-term of expansion (1381) in the sense of the pointweise convergence of 
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the space If g G T>(J\f) C in a similar way to equality (j 3 €3 H for (1 + n)t/i-order 

dual cumulant ( 1271) . n > 1, we derive 

lim i 2i+ +n (t, (y\x) 1; x) 5s _ n (y\x)^ = E M^ +n \z,X)g 8 _ n {Y\X)^ s = 

z c y\x, 

E^ E ^ ) (< 1 ,...,i fc ,x)p Ml (r\J0V'.. (4i) 

fe=l ii^...^i fc e{ji,...,jr s _„} 

Then, according to equalities (J7j) and (Hill , for group (J3"8"j) we obtain 
lim -((£/+ (t)g) s -g s ^ s = \imj(%i(t)g s -g s )ij s + 

E Tjr^i E Sg 7 a^ft (n*)i,*)^avw. = 

1 s 1 s 

AC^s + E^T E (b-r,\\ S •N£kUl>--->jk)9*-n[Y\{jl, ...,jn})A, 

I v • I Aj loft 

n=l k=n+l v y j l9 L...ft k =l 



where we used identity ( 12811 . 

Thus if g G £>(23 + ) C £ 7 (jF w ) in the sense of the *-weak convergence of the space 
£ 7 (^ r w ) we finally have 

w*-\im(±(U + (t)g-g)-<B + g)=0, 
where the generator Q3 + = ©^L 23^ of group of operators (138]) is given by ( |39ll . □ 

3.2 The existence and uniqueness theorem 

The following theorem holds for abstract initial- value problem (fl8|l - (fl9|l on the space 

Theorem 2. A solution of the initial-value problem to the quantum dual BBGKY hier- 
archy (Fl8|l - (1T9|1 is determined by the expansion 

1 

G s (t,Y) = J2 773^! E %t + n{t,(Y\X)i,X)G s -n(0,Y\X), s>l, (42) 

where the (1 + n)th-order dual cumulant %lf +n (t, (Y\X)i, X) is defined by (12711 
%t + nb( Y \ x h x ) ■= E (- 1 )' Phl (l p l " !) ! II *W*.*<)> 

P:{(Y\X)i,X}=U i X X,CP 

where J2p *' s the sum over all possible partitions P 0/ £/ie set {(Y\X)i, ji, . . . ,j s ~n} into 
|P| nonempty mutually disjoint subsets Xi C {(Y\X)i, X} . ForG(0) G r>(23 + ) C £ 7 (JF W ) 
it is a classical solution and for arbitrary initial data G(0) G £ 7 (JF W ) ii is a generalized 
solution. 
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Proof. According to Theorem 1, for the initial data G(0) G V( ( 3 + ) C £ 7 (.F W ), sequence 
( )42l) is a classical solution of initial- value problem ( |T8l) -(|T9l) in the sense of the *-weak 
convergence of the space £ 7 (^-^). 

Let us now show that in the general case G(0) G £ 7 (.F W ) expansions (142|) give a 
generalized solution of the initial-value problem to the quantum dual BBGKY hierarchy 
( TT5]) - (fT9|) . To this aim we consider the functional 

oo 1 

(f,G(t)) : ^-Tr ; s f s G s (t), (43) 

where / G £} a is a finite sequence of the degenerate trace class operators with infinitely 
times differentiable kernels and with compact support. According to estimate (j4"0l) this 
functional exists provided that a = 7 _1 > e (see Section 4). 
Using (14*21) . we can transform functional ()43l) as follows 

(/,G(f)) = (/,f/ + (t)G(0)) = ([/(*)/, G(0)). (44) 

In this equality the group £7 + (t) is defined by expression (1381) (Theorem 1) and U(t) is an 
adjoint mapping to the group U + (t) 

oo 

(C/(t)/) s (y) = V -:Tr s+lj ... jS+n 2t 1+n (t, F 1; A\F)/ s+n (A), (45) 

n=0 

where X = {1, . . . , s + n}, i.e. X\Y = {s + 1, . . . , s + n}. In expansion fT45]) the evolution 
operator 2l 1+n (t, Y 1; A\F) is the (1 + n)t/i-order cumulant of group of operators (IT5|) 

2i 1+n (t,r 1; = ^ (-i) |Phl (|P| - 1)! n 

P:{yi,x\y}=UiXi x z cp 

where is the sum over all possible partitions P of the set {Yx,X \ Y} = {Y\,s + 
1, . . . , s + n} into |P| nonempty mutually disjoint subsets Xi C {Y 1: X \ Y}. If / G Xi^, , 
series (l4"5j) converges, provided that a > e [14] and the functional (U(t)f,G(0)) exists. 

The one-parameter family of operators U(t) is differentiable with respect to t and for 
/ G £^0 an infinitesimal generator of group (j4"51) is defined by the following expression 
(see also ([22])) 

= -K(Y)f s (Y)+ 

E f:^!,..^^-^)^,...,^^)/^). (46) 

k=l ' i 1 ^z...^zi k =l n=l 

Since for bounded interaction potentials (TjJ, if / G £q , the operator *BU(t)f is a 
trace class operator, the operator %$U(t)fG(0) is also a trace class operator then the 
functional (<BE/(t)/,G(0)) exists. Moreover, it holds the equality: (<BU(t)f,G(0)) = 
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(C/(t)*8/,G(0)) and the following result 

gi|(i(^(t)-J)/,G(0))-(®/,G(0))| = 

oo 

gglE^i ^w(t(^(*)/ - /)/?.(«) - (»/) G.(0)) | < 

s=0 

oo 1 

ll G (°)lk(^) &S £ t-^wIt^)/ - /). - W)s\ = 

ll G (°)lk(^) &slli( l7 (*)/-/)-®/lk 7 _ 1 (^)= - 

Hence equality (jHJ) can be differentiated with respect to time and we get finally 

j t (f,G(t)) = (tf(f)®/,G(0)) = (*f,U + {t)G(0)) = 
(®/,G(f)), 

where the operator 03 is defined by (1461) . 

These equalities mean that the sequence of operators (j4~2l) for arbitrary G(Q) G £ 7 (^-^) 
is a generalized solution of the Cauchy problem to the quantum dual BBGKY hierarchy 

(HHD-dHD. □ 



4 The existence of the mean value observable 
functional 

As it was above mentioned, the functional of the mean value (fT5l) defines a duality be- 
tween marginal observables and marg inal states. If G(t) G £ 7 (.F W ) and F(0) G Z l a then, 
according to estimate (1401) . the functional 

oo ^ 

<G(t) |F(0)> = -1 Tr i,-,s G fl (t, 1, • • • , s)F,(0, l,...,s). (47) 

exists, provided that a = 7 _1 > e, and the following estimate holds 

|<G(t)|F(0)>|<e 2 (l- 7 e)- 1 ||G(0)|| £7( ^J|F(0)|| £; _ i( ^ ) . 

Thus, marginal density operators from the space describe finitely many quantum 
particles. Indeed, for such additive-type observable as number of particles, i.e. one- 
component sequence N(0) = (0, J, 0, . . . ) , according to definition ( }34"1) of dual cumulants, 
the expansion for solution ( 126]) get the following form (see also ( l29il ) 

s 

(N(t)) s (Y) =%t(t,l,...,s)J2 I = ^s,i, s > 1 

3=1 

and we have 

|<iV(t)|F(0)>| = iTnF^O,!)! < ||^(0)|| £ x_ i( ^ < oo. 
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We have above stated the properties of group ( l38i) defined on £ 7 (jF w ). In order to 
describe the evolution of infinitely many particles [1] the problem lies in the definition of 
functional (H7|) for operators from the suitable Banach spaces. Namely, marginal density 
operators have to belong to more general spaces than £* (JF W ). For example, it can be the 
space of sequences of bounded operators containing the equilibrium states [17]. In this 
case every term of expansions for the mean value functional (1471) contains the divergent 
traces [1], [2], [14] and the analysis of such a question for quantum systems remains an 
open problem. 

5 Conclusion 

The concept of dual cumulants (1541) of groups of operators §5§ for the Heisenberg equations 
([1]) or cumulants fj4"5]) of groups of operators ffl~3l) for the von Neumann equations ffTUl) 
forms the basis of the groups of operators for quantum evolution equations as well as 
the quantum dual BBGKY hierarchy and the BBGKY hierarchy for marginal density 
operators [14]. 

In the case of quantum systems of particles, obeying Fermi or Bose statistics, group 
of operators (1381) has different structures. The analysis of these cases will be given in a 
separate paper. 

On the space £ 7 (JF W ) one-parameter mapping (I38j) is not a strong continuous group. 
The group {U + (t)} t€ ^. of operators (1551) defined on the space £ 7 (JF W ) is dual to the strong 
continuous group {U(t)} t £R of operators (T4^|) for the BBGKY hierarchy defined on the 
space £^,(JF W ) and the fact that it is a Cg-group follows also from general theorems about 
properties of the dual semigroups [4,6]. 

We have constructed infinitesimal generator (1391) on the domain £>(*B + ) C £ 7 (.Fh) 
which is everywhere dense set for the *-weak topology of the space £ 7 (^>^). The question 
of how to define the domain of the definition T>( <! B + ) of generator (159"j) is an open problem 
[4,10]. 
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